Triatomic continuum resonances for large negative scattering lengths 
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We study triatomic systems in the regime of large negative scattering lengths which may be more 
favorable for the formation of condensed trimers in trapped ultracold monoatomic gases as the 
competition with the weakly bound dimers is absent. The manipulation of the scattering length can 
turn an excited weakly bound Efimov trimer into a continuum resonance. Its energy and width are 
described by universal scaling functions written in terms of the scattering length and the binding 
energy, B3, of the shallowest triatomic molecule. For a -1 < — 0.0297^/m Bi/% 2 the excited Efimov 
state turns into a continuum resonance. 

PACS numbers: 03.75.Fi, 36.40.-c, 34.10.+X, 21.45. +v 



O 

CO 



T3 
S3 
O 

> 

&\ 
(N 
(N 

m 
O 
-i— > 



I 

-a 
c 

o 
o 



X 



It is by now well established that shallow dimers are 
formed in trapped ultracold or condensed monoatomic 
gases, as it has been reported for 23 Na 87 Rb Q, and 
85 Rb 0. In the experiment in Ref. Q, an atom-molecule 
coherence in the Bose-Einstein condensate was also ob- 
served. The measured oscillation frequency of the quan- 
tum superposition of 85 Rb dimers and atoms in the con- 
densate was in agreement with the shallow 85 Rb2 binding 
energy over a wide range of values near a Feshbach res- 
onance. However, the formation of triatomic molecules 
has not yet been observed. 

Recently, the energy of the shallowest bound tri- 
atomic molecule in trapped ultracold and condensed 
monoatomic gases was predicted using the value of the 
measured recombination rate into a shallow dimer; the 
energies ranged from 7.75 mK down to 0.24 nK near a 
Feshbach resonance Q. Condensed triatomic molecules 
coexisting with dimers and atoms near a Feshbach res- 
onance would present an oscillatory dependence of ob- 
servables on the trimer binding energy [5j. For negative 
scattering lengths and zero energy, the recombination 
rate into deep diatomic molecular states shows a resonant 
peak at values of a for which the trimer Efimov state 
hits the three-body continuum threshold • If the mag- 
nitude of the large negative scattering length is decreased 
after an Efimov state hits the continuum, it turns into a 
three-body resonance, as we will show. Therefore, at non 
zero energies, or temperatures, the resonant peak of the 
recombination rate would in principle appears when the 
resonance energy matches the energy of the continuum 
triatomic system. 

Furthermore, the manipulation of the scattering length 
in monoatomic condensates in the regime of large and 
negative values offers an interesting possibility. No shal- 
low dimers can exist in this case. Then, the formation of 
triatomic molecules may be more favorable as the main 
competitors are absent. If a coherent quantum super- 
position of atoms and trimers or resonances appears, it 
will present an oscillatory frequency corresponding only 
to the shallow three-body bound or resonant state. As 
we will show, this frequency scales with the scattering 



length and with a triatomic physical scale in the form 
of a universal function. Here, the binding energy of the 
shallowest trimer state is chosen as the three-body phys- 
ical scale. 

In 1970, Efimov predicted that infinitely many weakly 
bound three-boson states appear when the s— wave two- 
boson scattering lengths, a, goes to the limit of a — > 
±00 |fj. For large scattering lengths, an attractive long- 
ranged effective interaction binds the three-particle sys- 
tem in a range of about \a\ (See also 0)- A new bound 
state appears for every increase of \a\ to ~ 23|a|. For 
positive a, at the threshold of the new state, the trimer 

binding energy is 6.9B 2 s mini (£ 2 is the dimer bind- 
ing energy). For a < 0, a new state becomes bound when 
the trimer has an energy of - 1100.B2 @, where now Bo 
is the energy of the virtual dimer state. Due to the large 
size of the system, the threshold conditions for the ex- 
istence of excited Efimov trimers are universal 0, Il2| . 
i.e., independent of the detailed potential shape, and ex- 
hibit a scaling form for large values of a/r [EM El (the 
interaction range is ro). The three-body system heals 
through regions that are outside the potential action, 
where the wave function is essentially a solution of the 
free Schrodinger equation, and therefore the properties 
of the system are defined by few physical scales. 

The dimer and trimer binding energies are the only 
physical scales that survive in the limit of (a/rg) — ► ±00 
(scaling limit), which essentially relates the Thomas col- 
lapse of the trimer state for r$ — > 01 to the Efimov ef- 
fect (\a\ -> 00) [l|. In the scaling limit, the three-boson 
observables are functions of the shallowest trimer bind- 
ing energy (reference three-body energy) and B 2 . These 
functions approach universal curves |y, |9( . The collapse 
of the three-boson system in the limit of a zero-range 
force makes the trimer energy the three-body scale of 
the system beyond the two-body energy [l5j |. 

For large scattering lengths, an excited Efimov state 
turns into a virtual state when a > is decreased |16). 
The threshold moves faster than the energy of the excited 
state. Therefore, with the increase of a > 0, states pump 
out from the second sheet of energy to become bound 



2 



states [10|. If a < is decreased in magnitude, the trimer 
bound states dive into the continuum 9J. It is our aim 
here to evaluate the scaling properties of the energy and 
width (for the decay into the three-body continuum) of 
the resonance born from an Efimov state, when a large 
a < is varied. 

The 4 He excited trimer state calculated also with real- 
istic models offers a good example of an Efimov state and 
its universal scaling properties with the shallow dimer 
and trimer binding energies 0, ITU . These molecules are 
special due to the large spatial size which spreads out 
much beyond the range of the potential [HEl]. The 4 He- 
4 He root-mean-square distance in 4 He3 for the ground 
and excited states 0, are of the order of 5 to 10 A 
and of about 50 to 90 A, respectively. The product of 
the mean-square interatom distance with the separation 
energy of one atom from the trimer in units of % = m = 1 
(m is the atom m ass) is about the 1 in the ground and 
the excited states [171 llSf , while the ratio of the binding 
energies /Bjp » 500 |(| (-Bg ' is the binding energy 
of the N th trimer state). 

In the present work we calculate the three-boson res- 
onance energy and width for the decay into continuum 
states for large and negative scattering lengths. In this 
case it is justifiable to use a Dirac-<5 potential. We 
solve subtracted homogeneous equations defined within 
a renormalization scheme applied to a three-body sys- 
tem interacting with s— wave zero-range pairwise poten- 
tials [HEi. We present the results for the energies and 
widths in the form of a universal scaling function 0, 
which gives the trajectory of the three-body energy in 
the complex plane as a function of the shallow two-body 
virtual state energy. We show that a resonance becomes 
an excited trimer Efimov state when the large a < is 
decreased. We go beyond Ref. where it was found 
the dependence of the virtual three-boson state energy 
with a > born from an Efimov state which entered in 
the second energy sheet. 

One can fix the energy of one three-body bound-state 
(the three-body physical scale), and the two-body scat- 
tering length and get other observables. All the detailed 
information about the short-range force, beyond the low- 
energy two-body observables, is retained in only one 
three-body physical information in the limit of zero-range 
interaction. The existence of a three-body scale implies in 
the low energy universality found in three-body systems, 
or correlations between three-body observables [15J, |20( . 
In the scaling limit 0, E| j one nas 

O (E, B 3 , B 2 ) (B 3 )-« = T ^E/B 3 ,±y/B 2 /Bs) , (1) 

where O is a general observable of the three-body system 
at an energy E, with dimension of energy to the power r\. 
This equation means that any observable of the system 
can be represented by a function that depends only on 
one three- and one two-body scales. The three-body scale 
is brought by the reference energy B% (binding energy) of 
the shallowest trimer state. The ± sign denotes positive 



or negative scattering lengths. In the case of the energies 
of a resonance, of an excited or virtual trimer states, 
instead of Eq. , the scaling function is written as: 

E 3 = B 3 £ (±y/B 2 /B 3 ) . (2) 

(Throughout this paper, we use units such that h = m = 
1, where m is the mass of the atom.) 

After partial wave projection, the s— wave subtracted 
integral equation for three identical bosons is given 

by m 

X {q)=^r{£) J dq'q' 2 I dy X (q') (3) 

\E 3 - q 2 - q' 2 - qq'y -(-i 2 - q 2 - q' 2 - qq'yj ' 

where £ = E 3 — |g 2 , fi is the subtraction point and r is 
given by 

r-\0 = -2**y/Bl-*xt , (4) 

Jo t-P 

here B 2 is the dimer virtual state energy. (For positive 
values of the real part of E 3 , Eqs. 10 and (0J are analyt- 
ically extended to the second energy sheet, as we discuss 
below.) 

Throughout this work we perform calculations only 
considering a < for a virtual dimer state. We use a 
contour deformation method to calculate the resonance 
energy and width [2l|. The homogeneous equation, 
is analytically continued to the second sheet of energy, 
by making q (q 1 ) — > q e~ %B (q'e^ 8 ) with < 6 < 7r/4. 
For large enough 8, the solution of Eq.© in the complex 
energy plane is found for tan(26>) > —Im(E 3 )/Re(E 3 ). 

In the limit of a virtual dimer energy tending to zero an 
infinite number of Efimov states appears from the solu- 
tion of Eq. © for negative E 3 (= -B 3 ). The N th Efimov 
state has binding energy given by B$ N ^ with N = indi- 
cating the ground state obtained from Eq. © . When the 
N th Efimov state turns into a resonance its complex en- 
ergy is denoted by -E^^ . In figure , it is shown the real 
part of the complex energies, in units of fi — 1, for the 
first three states obtained by solving numerically Eq. ||3J| . 
In the figure, the values of [Re(E 3 )] 1 ' 2 for the first three 
resonances are described by the positive part of the plot. 
The curves represent -Eg ^ (solid line), E^ (dashed line) 

(2) 

and £"g ; (dotted line), i.e., the complex energies of the 
first three states of Eq. ©. Note that the subtraction 
present in Eq. J3J) regularizes it, and the Thomas collapse 
is avoided. When B 2 is decreased the resonance turns 
into a bound state and the negative part of the plot give 
the energies of the ground (-Bg ), first (-Bg ) and second 

(Bg ) bound trimer states. One realizes that the form 
of the curves are very similar which indicates that the 
function £ of Eq. does not depend on the state, which 
will be confirmed later on, in the scaling plot of figure ©• 
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FIG. 2: Imaginary part of the three-body resonance energy 
as a function of the two-body virtual state energy in units of 
fi = 1. Curves labelled as in fig.l . 

picture of the route of an Efimov state when a is var- 
ied passing through a Feshbach resonance. If one begins 
with positive a and increases it, a virtual trimer state 
becomes bound, then crossing a Feshbach resonance and 
decreasing the magnitude of the large negative value of 
a, the trimer bound state becomes unbound and turns 
into a resonance. The results for the trimer bound state 



FIG. 1: Real part of the three-body resonance energy as a 
function of the two-body virtual state energy in units of fj, = 1. 
The negative values represents the results for square-root of 
the energy of the N th bound trimer state, [B^^] 1 ^ 2 , which 
turns into a resonance when \B2\ is increased. The positive 
values of [Re(E 3 N ' > ] 1 ^ 2 corresponds to the resonance. The 
results are given for TV = (solid line), N = 1 (dashed line) 
and N = 2 (dotted line). 



The values of v ' B-2 in units of /j, = 1 at which the Efimov- 
state becomes unbound are 0.066, 0.0028 and 0.00013 for 
the ground, first and second states, respectively. The 
ratios (0.066/0.0028) 2 « (0.0028/0.00013) 2 - 500 are 
practically independent of the state, in agreement with 
the scaling limit implied by Eq. (J2J. It is curious that 
the real part of the resonance energy tends to zero for B2 
large enough, while the width is nonzero. 

In figure @ the results for the imaginary part of the 
resonance complex energy are shown as a function of B2 
in units of fi = 1. The solid, dashed and dotted lines are, 
respectively, the correspondent imaginary part of E^\ 
E^ and E$ . The threshold values of B2 for which 



the resonant state becomes bound are clearly seen in the 
figure and the resonance width, rl) = 2| Im[E^ N ']\, in- 
creases with the virtual dimer energy. The results shown 
in figures Q and @ give the complete trajectory of a 
three-body bound Efimov state when the two-atom inter- 
action with a < is changed, as can be done for trapped 
atoms near a Feshbach resonance. These results extend 
the findings of Ref. 10]. Now, we can give a complete 
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FIG. 3: Ratio of the real part of E { 3 N) and B. 



,(JV-1) 



as a func- 



(JV-l) 



The solid line is the results for N=l and 



tion of B 2 /B. 
the dashed line for N=2. 

or resonance energies can be expressed in the form of a 
universal scaling function, Eq. J5J), that depends only on 
the ratio of the two- and the reference three-body ener- 
gies. From the calculations given in figures and 
we construct the scaling plots of figures © and Q}, re- 
spectively. For this purpose, we use the first and second 
Efimov states (bound or resonant) and we numerically 
obtain 
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using N equal to 1 and 2. (In Eq. JSJ the states are in 
fact indexed, which was not explicitly shown in Eq. J2J))- 
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FIG. 4: Ratio of the imaginary part of E$ and Bg 1 as 
a function of B2/-B3 -1 . Curves labelled as in fig. 3 . 

We show the results for the real and imaginary part of 
the resonance energy in a form of scaling plots in figures 
© and Q), respectively. In these figures the solid line 
corresponds to N = 1 and the dotted line to N = 2. 
In figure JSJl, the negative values of [Re(E 3 )/ -B3] 1 / 2 give 
the trimer bound state results, while the positive val- 
ues come from the resonant state. It is clear that a 
universal curve for the scaling function £(—^/B 2 /B 3 ) is 
approached as the results for N = 1 and 2 practically 
overlap. The imaginary part of the resonance energy 
in figure (J3J appears for the dimer virtual state energy 
B2 > 0.000882 B3, where the excited Efimov bound state 



turns into th e resona nce. Our results extends the scaling 
function £(-J B^/B^) (bound dimer states) obtained pre- 
viously in |3| and [13 to negative values of the argument 
and for the trimer resonance region. 

In conclusion, we have obtained the trajectory of an 
Efimov state when the two-atom interaction crosses a 
Feshbach resonance, in a form of a scaling function. The 
value of the resonance energy and width depends only 
on the dimer bound or virtual energy and the shallowest 
trimer binding energy. The route of an Efimov state when 
a is varied passing through a Feshbach resonance can be 
summarized as follows: beginning from a large positive a 
and increasing it further, a virtual trimer state becomes 
bound and then after crossing a Feshbach resonance and 
decreasing the magnitude of the large negative value of a, 
the trimer bound state becomes u nbound a nd turns into 

a resonance, for a -1 < — 0.0297v/m B3/F1 2 , where B3 is 
the binding energy of the shallowest trimer state. In this 
respect, the study of trapped Bose-Einstein condensates 
near a Feshbach resonance can be fruitful not only to 
reveal the properties of an interacting coherent quantum 
state under extreme conditions, but could also give new 
insights in the few-body physics underlying the curious 
phenomenon of Efimov resonant states. 

We would like to thank Prof. Y. Koike from Hosei Uni- 
versity for a helpful discussion. We thank the Brazilian 
agencies FAPESP (Fundacao de Amparo a Pesquisa do 
Estado de Sao Paulo) and CNPq (Conselho Nacional de 
Desenvolvimento Cientffico e Tecnologico) for financial 
support. 



[f] C. Mckenzie et al, Phys. Rev. Lett. 88, f 20403 (2002). 

[2] R.H. Wynar, R.S. Freeland, D.J. Han, C. Ryu, and D.J. 
Heinzen, Science 287, f0f6 (2000). 

[3] N.R. Claussen, E.A. Donley, S.T. Thompson, and C.E. 
Wieman, Phys. Rev. Lett. 89, 01040f (2002); E.A. Don- 
ley, N.R. Claussen, S.T. Thompson, and C.E. Wieman, 
Nature (London) 417, 529 (2002). 

[4] M.T. Yamashita, T. Frederico, A. Delfino, L. Tomio, 
Phys. Rev. A68, 033406 (2003). 

[5] E. Braaten, H.-W. Hammer, M. Kusunoki, Phys. Rev. 
Lett. 90, 170402 (2003) 

[6] V. Efimov, Phys. Lett. B 33, 563 (1970); Nucl. Phys. 
A362, 45 (1981); V. Efimov, Comm. Nucl. Part. Phys. 
19, 271 (1990) and references therein. 

[7] E. Braaten, H.-W. Hammer, Phys. Rev. Lett. 87, 160407 
(2001). 

[8] E. Nielsen, D. V. Fedorov, A. S. Jensen, E. Garrido, Phys. 

Rep. 347, 373 (2001). 
[9] A. E.A. Amorim, T. Frederico and L. Tomio, Phys. Rev. 
C56, R2378 (1997); T. Frederico, L. Tomio, A. Delfino 
and A. E.A. Amorim, Phys. Rev. A60, R9 (1999). 
[10] M.T. Yamashita, T. Frederico, A. Delfino, L. Tomio, 

Phys. Rev. A66, 052702 (2002). 
[11] E. Braaten, H.-W. Hammer and M. Kuzunoki, Phys. 



Rev. A67, 022505 (2003). 
[12] A. Delfino, T. Frederico, L. Tomio, Few-Body Syst. 28, 

259 (2000); J. Chem. Phys. 113, 7874 (2000). 
[13] L.H. Thomas, Phys. Rev. 47, 903 (1935). 
[14] S.K. Adhikari, A. Delfino, T. Frederico, I.D. Goldman, 

L. Tomio, Phys. Rev. A37, 3666 (1988). 
[15] S.K. Adhikari, T. Frederico and I.D. Goldman, Phys. 

Rev. Lett. 74, 487 (1995); S.K. Adhikari and T. Fred- 
erico, Phys. Rev. Lett. 74, 4572(1995). 
[16] S. K. Adhikari, A. C. Fonseca and L. Tomio, Phys. Rev. 

C26, 77 (1982); S. K. Adhikari, L. Tomio, Phys. Rev. 

C26, 83 (1982). 
[17] P. Barletta and A. Kievsky, Phys. Rev. A64, 042514 

(2001). 

[18] V. Roudnev, S. Yakovlev, Chem. Phys. Lett. 328, 97 
(2000). 

[19] M.T. Yamashita, R.S. Marques de Carvalho, L. Tomio 
and T. Frederico, Phys. Rev. A68, 012506 (2003). 

[20] T. Frederico, I.D. Goldman, Phys. Rev. C36 1661 (1987). 

[21] J. Aguilar and J.M. Combes, Commun. Math. Phys. 
22, 269 (1971); E. Balslev and J.M. Combes, Commun. 
Math. Phys. 22, 280 (1971). 



